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Abstract 

We calculate analytically the distributions of "level curvatures" (LC) (the 
second derivatives of eigenvalues with respect to a magnetic flux) for a particle 
moving 

in a white-noise random potential. We find that the Zakrzewski-Delande 

conjecture is still valid even if the lowest weak localization corrections 

are taken into account. The ratio of mean level curvature modulus to mean 

dissipative conductance is proved to be universal and equal to 2n in agreement 

with available numerical data. 
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Nowadays it is considered to be a well established fact that spectral statistics of generic 
chaotic or disordered systems is adequately described by that typical for eigenvalues of 
large Gaussian random matrices (RM) |I|]. The range of applicability of the "Gaussian 
universality" to real disordered systems was considered in the pioneering papers by Efetov 
and Altshuler and Shklovski . It was demonstrated that statistical properties of the energy 
levels of a quantum particle moving in a static random potential follow RM predictions as 
long as effects of Anderson localization on the particle diffusion are neglegible. Another 
important result is due to Berry who demonstrated how the same universality may arise in 
globally chaotic ballistic systems ( "quantum billiards") 0]. 

Quite recently an interesting new developement in the study of spectra of disordered 
systems and their chaotic counterparts has been made. The problem is to study the so- 
called "level response statistics", i.e. to provide a statistical description of sensitivity of the 
energy levels to external perturbations of different types. This issue attracted a great deal of 
research interest, both analytically and numerically [^]- 0, [|11|]- The most frequently 



studied characteristics are first and second derivatives of the energy levels En{a) with respect 
to a tunable parameter a characterizing the strength of perturbation. Physically the role 
of such a parameter can be played by, e.g. an external magnetic field, the strength of a 
scattering potential for disordered metal, a form of confining potential for quantum billiards, 
or any other appropriate parameter on which the system Hamiltonian is dependent. The 
first derivatives f„ = dEn/da are frequently called the "level velocities" (LV) (or "level 
currents"), the second ones = d'^En/da^ are known as "level curvatures". 

In a series of papers by the MIT group see also [0], it was found that for a generic 
chaotic system whose unperturbed spectrum is well described by the universal RM statistics 
the set of "level velocities" fn(a) is characterized (after appropriate rescaling) by a universal 
correlation function {vn{a)Vn{o.')) whose form is again dependent only on the symmetry of 
the unperturbed Hamiltonian and that of the perturbation. 

The range of applicability of these results to real systems is the same as before: they 
are valid for systems with completely "ergodic" chaotic eigenf unctions covering randomly. 



but uniformly all the available phase space and showing no specific internal structure. The 
universality is believed to be independent of "whether the chaos originates in mesoscopic 
disorder or deterministic instability of the classical trajectories." M. The effects of eigenfunc- 



tion localization- either due to scarring |T^ or due to disorder-induced quantum interference 
(the Anderson localization) - result in substantial modifications of the LV characteristics, 



see 



Much interest was concentrated on the "level curvature" (LC) distribution V{K). Gas- 
pard and coworkers Jill discovered the universal asymptotic behaviour V{K) oc in 



the large curvature limit K —>■ oo, the parameter /3 = 1, 2 or 4 depending on the symmetry 
universality class. This behaviour is a direct consequence of the so-called "level repulsion" 
for unperturbed systems. Namely, the probability V{s) for two unperturbed energy levels 
to be separated by the spacing s much smaller than the mean level spacing A vanishes 
as V{s) oc s^, and this fact can be shown to result in the abovementioned large-curvature 
behaviour W^. The analytical form of the whole distribution function V{K) was guessed by 



Zakrzewski and Delande Wl\ on the basis of numerical results and is given by the following 



expression: 



p(fc) = C(/5)(l + A;2)-{i+/^/2) (1) 

where Cp is a normalization constant, /? = 1, 2, 4 depending on presence or absense of time- 
reversal invariance and Kramers degeneracy and k is the dimensionless level curvature. We 
will refer to this expression as to the ZD conjecture. 

Very recently, von Oppen succeeded in demonstrating that eq.(|lD is indeed exact 



for the ensembles of large Gaussian Hermitian {(3 = 2) and real symmetric matrices. The 
validity of ZD conjecture for all three classes of large Gaussian matrices was proved in a 



different way by the present authors [|2^. It is natural to suppose, that the status of the 
distribution eq.(|I|) within the theory of disordered and chaotic systems is the same as before: 
it is valid for systems with completely ergodic extended eigenstates. This was indeed found 
to be the case in a series of interesting numerical experiments ||T9|j20[| on quasi-lD as well 



as 3D periodic random tight-binding models subject to the influence of Aharonov-Bohm 
magnetic flux, which acts as a time-reversal symmetry breaking parameter. In particular, 
the distribution eq.(|l]) was found to persist up to the Anderson localization transition |I9| . 



Unfortunately, the methods used in 1T^J2^ for deriving the form of LC distributions for 



the Gaussian ensembles were havily based upon the explicit form of the joint probability 
density of eigenvalues known for all these ensembles f^, which is of course an immense 
simpliflcation. One therefore has to invent a different technique in order to be able to treat 
the curvature distribution analytically under more realistic assumptions. 

It turns out to be possible to flnd such a technique for the case of time-reversal invariant 
systems subject to a time-reversal symmetry breaking perturbation. This case seems to be 
one of the most interesting from the physical point of view as well as relevant experimentally. 
As a physical realisation one can imagine a disordered mesoscopic sample (e.g. cylinder or 
ring) pierced by magnetic flux 0. For such a system the "typical" level curvature is expected 
to be related to the dimensionless conductance gc of the sample due to the famous Thouless 
formula: ■^d'^En/d(j)'^\^=o ~ Ec/ A = gc- Initially suggested by Thouless this relation 
attracted renewed interest recently. Its meaning was reconsidered in a broader context by 
Akkermans and Montambaux |2^, see also a quite detailed discussion in [^,^. All these 



facts make the consideration of such systems to be of special interest. 

A speciflc feature allowing us to treat the particular case of weak time-reversal symmetry 
breaking perturbation analytically is vanishing of the flrst derivatives dEn/ d(j)\^=Q on reasons 
of symmetry. This fact allows one to represent the curvature distribution in terms of a 
product of advanced and retarded Green functions, the average over the disorder being 
performed nonperturbatively with help of Efetov's supersymmetry approach [Q] . The method 
provides a unique possibility to derive the level curvature distribution starting from a genuine 
microscopic Hamiltonian of a quantum particle experiencing elastic scattering: 

n{ct>) = ^{p--^A{ct>))\u{r) (2) 

with f/(r) being a white noise random potential and A(</)) standing for the vector potential 



corresponding to the magnetic flux 0. 

This fact allows one to try to take into account the weak localization effects due to the 
flnite ratio of the Thouless energy Ec to the mean level spacing A. The general way of doing 



this was developed recently by Kravtsov and Mirlin Exploiting this method one can 
flnd that the curvature distribution eq.(|I|) preserves its form to the flrst order in A/ Ec <ti I, 
the width {\K\) being renormalized. 

To begin with, we introduce the resolvent operator G^{a; e) = [E — T[{(f)) ± where 
a = 0/00, with 00 = 27ic/e being the flux quanta. Let us now consider the following 
correlation function: 

/C(m) = 



lim,^o eTrG+(a = 2^e/u; e)TrG-(a = 0; e) = (3) 
lim - 

e^O Z^„,m=l [E-E„{a=0)-ie][E-Em{a=2y/7/Z)+ie] 

Here we used the expression for the trace of the resolvent TrG in terms of eigenvalues of the 
Hamiltonian H{4>). 

It is important to note that the limiting procedure e ^ in eq. @ is performed prior to 
the thermodynamic limit V — * oo, where V is the system volume. Therefore, the parameter 
e (that plays a role of effective "level broadening" necessary to regularize the resolvent 
operator) can be considered as small in comparison with the mean level spacing A oc 1/^. 
It was already mentioned that the probability of having two levels at a distance s <^ A is 
vanishingly small due to level repulsion. Taking this fact into account one flnds that the only 
terms that survive the limiting procedure e — > are those with coinciding indices m = n. 
Remembering dEn/d(f)\(f,=o, one obtains: 

>^(-)--i:-r^i(E-E„y. /f» = ^l„o (4) 

Let us perform formally the averaging over the disorder and introduce the function V{K) = 
A J2n=i ^(-^ ~ Kn)S{E — Enf^. This function has the meaning of distribution of curvatures 
for levels in a narrow spectral window around the energy E . Then one obtains the following 
relation: 



-{K{u))=u /"dwe-"" 

TT JO 



dKViK)exp {i vK) 



(5) 



Let us stress that it is simultaneous vanishing of all the first derivatives f„ = dEn/dcf) 
at = that made it possible to scale 0/00 oc e^/^ producing the finite limiting expression 
eq.(^) when e ^ in eq.(^. In contrast, when f „ 7^ one can scale oc e and get 
the expression eq.(H), but with the level velocity distribution V{v) substituted for the level 
curvature distribution V{K). This fact was already used in p!2|j8|] in order to calculate 
analytically LV distributions for different systems. 

In order to evaluate the average correlation function /C(m) one can use Efetov's supersym- 
metry approach 0]. Exploiting the weak disorder parameter kpl ^ 1, with kp denoting the 
Fermi wavenumber and I denoting the mean free path due to elastic scattering, the problem 
can be mapped onto the so-called nonlinear graded a- model. Rather detailed exposition 



of the mapping can be found, e.g. in p5[. The final expressions appeared frequently in 



literature, most recently in p^, and can be used for our needs after a slight modification 
resulting from the fact that in our case the magnetic flux enters only one of two Green 



functions rather than both as in [25,24|. As the result, one finds: 



(/C(n)) = ^lim,^oi./2^/i(Q)x 
/ f Str [g(r)A;i±A] / ^str [Q{r)k^-^] exp ~S,{Q) 

S.{Q) = 



(6) 



ttD r dr 



/fstr vg 



g,r 



7re r dr 

2A J y 



/fStrQA 



where D is the classical diffusion constant D oc vpl due to random scattering , is the 
vector potential corresponding to the magnetic flux through the sample 0/0o 



2ye/u in 

the gauge ensuring that A is tangential to the surface of sample; e.g. ^A = ^-^^t ^ot 
the ring geometry. The notation [ , ] stands for the matrix commutator and Str stands for 
the graded trace. The position dependent 8x8 supermatrices Q{r) satisfy the constraint 
Strg2=i and can be parametrized as follows: Q = T ^AT, where T belongs to a graded 
coset space U0SP(2, 2/4)/UOSP(2/2) x UOSP(2/,2). Other 8x8 supermatrices entering 
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the eq.(^) are diagonal: 

A = diag(/2, -h, -h) r = diag(o-, a, 0, 0, 0, 0) 
k = diag(J2, -/2, h, -h) = diag(l, -1); h = diag(l, 1) 
Let us now recall that there are two relevant energy scales for a quantum particle diffusing 
in a closed disordered sample of the size L > I 

: the mean level spacing A and the Thouless energy Ec = D/L"^. For A ^ i?c it is 
known ^ that 

the main contribution to the functional integral, eq.(|^) comes from spatially uniform 
configurations ( a so-called "zero mode"): Q{r) = Qq = T^^AT. The condition gc = 
Ec/A ^ 1 determines the regime where statistical characteristics both of eigenfunctions 
and energy levels of the Hamiltonian eq.(@) are adequately described by classical 

RM theory 0J^,|1^. The corrections to RM results have a form of regular expansion in 
the small parameter g^^. A systematic way to construct such an expansion was originally 
suggested in and discussed in more details in [|13|. The general procedure is as follows: 



one decomposes the matrix Q{r) as Q{r) = TQ'^Q{r)TQ and uses the fact that supermatri- 
ces Q{r) only weakly fluctuate around the value Q = A as long as gc ^ 1- As the result, 
the integration over Q{r) in eq.(|^) can be performed perturbatively. After that one ob- 
tains a renormalized expression for the "zero mode effective action" S^{Qo) governing the 
integral over spatially uniform configuration Qq. This last integration has to be performed 
nonperturbatively. 

Performing the "perturbative" part of the evaluation of the functional 
integral to the first order in the small parameter g~^ and keeping only terms relevant 
at e — i> one obtains for the quantity {IC{u)) the same expression eq,(|^) with the following 
replacement Q{r) Qq, Se{Q) Se{Qo), where 

-^Str(QoA) - ^bEc (l - / f n(r, r)) Str(gof Qot) 
Here we introduced a sample- dependent geometrical factor b = (^^^^ L'^ J y (f ^e) 5 
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[p^] , equal for the ring geometry to 6 = 47r^. In the expression eq.(|^) the quantity 
n(r, r') stands for the so called diffusion propagator whose Fourier transform is equal to 
^ ^ :, momentum q going over all nonzero values allowed from the quantization condition 



in a given sample geometry, v being the mean density of states at Fermi level. 

In order to evaluate the integral over Qq nonperturbatively one has to use the so-called 
Efetov parametrisation of the corresponding coset space. A useful discussion of necessary 



technical details can be found, in particular, in the paper ||2^. Below we present only major 
steps of the evaluation procedure relegating most of details to more extended publication 

After standard manipulations, the integration over the graded coset space is reduced to 
performing the following three-fold integral: 

^ {K{u)) = 

"l - f (A? + Xl + X^- 2AA1A2 - 1)] X 
exp [-27re(AiA2 - A) - n^ (A? + A^ - A^ - 1)" 

where we introduced the notations u = u (^^[1 — / ^n(r, r)]^ , e = 

A close inspection of this expression makes it clear that the only contribution nonvan- 
ishing in the limit e — is that coming from the region Ai ~ A2 oc e^^/^. Picking up such a 
contribution one gets after some algebra: 

- {K{u)) = 

2 A (9) 

2^ Jo "^^^ "Jo ' (^y+y-i-2X) 

Two last integrations in this expression can be performed analytically. Comparing the result 
with the relation eq. (|^) one obtains the following expression for the Fourier transform of the 
LC distribution: 

V{K) expivKdK = ^vKi{-iv), (10) 
-00 2, 

7 = 2bE,(^l- JyUir,r)j (11) 



where Ki{z) is the Macdonald function. Inverting this Fourier transform one finally arrives 
at the desired "level curvature" distribution: 

1 7^ 

-pfK) = 1 — (12) 

Thus, we proved that the level curvature distribution for weakly disordered Aharonov- 
Bohm systems follows 



exactly the form suggested by Zakrzewski and Delande ||Tl| even beyond the domain of 
validity of random matrix theory, i.e. when the lowest weak localization corrections are 
taken into account. This conclusion is in good coordination with available data from recent 



papers ||T9|J20[| where numerical simulations of lattice analogues of disordered Aharonov- 
Bohm systems were performed. Having in mind a connection betwen asymptotic behaviour 
of LC distribution and degree of level repulsion, it is worth mentioning that our result is 



also in agreement with the recent paper p3[ where the two-level correlator was shown to be 



affected by weak localization corrections only at the second order with respect to g~^. 

Another interesting point is that our exact results when combined with the earlier findings 
I^J^ allow one to prove the universality of the ratio of the LC mean scaled modulus (| -ft' | 
)/A = 7/A to the sample dissipative conductance P J2(]| , |2^ defined as C(0) = X2^(^f) 
where the average goes both over the flux value and over the disorder. The calculation done 
in produces the value C(0) = bEc/iTcA) for any Aharonov-Bohm diffusive system when 
(7c ^ 1- Comparing this result with eq.(ll) one finds ^J^.^^ l3c>i= 27r in a good agreement 



with the value 6.7 ± 10% found in numerical simulations pO . 

Finally, let us mention that our method allows one to extend the validity of the ZD 
conjecture to large symmetric RM with independent arbitrarily 

distributed entries having finite variance. In particular, the most interesting ensemble of 



large sparse RM |28| can be treated analytically in this way. These results will be published 



elsewhere 127 
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